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In  active  sonar  applications,  a  quadrature  receiver  can  be  employed 
to  detect  the  presence  of  an  underwater  reflector.  This  well-known 
receiver  performs  optimally  when  the  input  waveform  is  a  coherent  signal 
(i.e.,  constant  but  unknown  amplitude  and  phase)  embedded  in  white, 
zero-mean  noise.  Unfortunately,  when  an  attempt  is  made  to  detect  a 
distributed  highlight  reflector  in  a  multipath  environment  by  using  a 
narrowband  transmit  signal,  the  various  multipath  returns  may  overlap,  — 
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producing  a  non-resolvable  multipath  situation  and  causing  the  amplitude 
and  phase  of  the  signal  to  become  random  functions  of  time.  Under  these 
circumstances,  the  detection  performance  of  the  quadrature  receiver  is 
degraded  compared  to  its  performance  with  a  coherent  signal. 

This  study  theoretically  measures  the  degradation  in  the  receiver's 
output  signal-to-noise  ratio  (SNR)  when  the  input  waveform  consists  of 
a  narrowband  Gaussian  signal  embedded  in  white  noise.  The  output  SNR's 
dependence  on  the  correlation  characteristics  of  the  input  signal  is 
clearly  shown. 
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PERFORMANCE  DEGRADATION  IN  A  QUADRATURE  RECEIVER 
FOR  CW  SIGNALS  CORRUPTED  BY  MULTIPATH 


INTRODUCTION 

In  active  sonar  applications ,  a  quadrature  receiver  can  be  used  to 
detect  the  presence  of  an  underwater  reflector.  This  receiver's 
performance  is  optimum  when  the  amplitude  and  phase  of  the  reflector's 
signal  are  constant  but  unknown  (i.e.,  coherent  signal)  and  the  back¬ 
ground  noise  is  white  with  zero-mean.  As  shown  in  Fig.  1,  this  receiver 
correlates  the  received  waveform  (either  noise-only  or  signal-plus-noise) 
with  sine  and  cosine  waves  having  fixed  parameters;  the  respective 
correlator  outputs  are  then  squared  and  summed.  A  detection  is  made  if 
the  output,  Q(T) ,  exceeds  a  predetermined  threshold.  Figure  2  shows  an 
equivalent  form  of  the  quadrature  receiver,  the  familiar  matched  filter- 
envelope  detector  combination  referred  to  as  the  incoherent  matched 
filter  C13 •  In  many  sonar  applications  a  bank  of  matched  filter-envelope 
detectors,  each  centered  at  a  different  frequency,  is  used  to  detect  an 
underwater  reflector  and  estimate  the  reflector's  Doppler  frequency  and 
range . 

One  problem  in  active  sonar  systems  with  narrowband  transmit  signals 
is  the  combination  at  the  receiver  of  direct  path  and  surface  path  returns 
reflected  from  a  distributed  highlight  reflector.  When  these  multipath 
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s(t),  RECEIVED  SIGNAL 
n(t),  BACKGROUND  NOISE 
T,  INTEGRATION  TIME 

Figure  1.  Quadrature  Receiver 


r(t)«s(t)+n(t) 


h(t)=sin  2vf(T-t) 
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s(t),  RECEIVED  SIGNAL 
n(t),  BACKGROUND  NOISE 
T,  SIGNAL  DURATION 


Figure  2.  Incoherent  Matched  Filter 


returns  overlap  (i.e.,  non-resolvable  multipath),  the  composite  echo  has  a 
random  time-varying  amplitude  and  phase  over  the  signal  duration.  The 
mismatch  that  exists  between  the  random  time-varying  parameters  of  the 
echo  and  the  fixed  parameters  of  the  sine  and  cosine  multipliers  in  the 
quadrature  receiver  causes  a  degradation  in  detection  performance.  This 
degradation  has  been  observed  and  quantified  empirically  but  is  not  well 
understood  theoretically. 

The  problem  of  performance  degradation  in  matched  filter  or  replica- 
correlator  processors  has  been  dealt  with  previously  by  several 
researchers.  Costas  [2]  examined  the  detection  performance  losses  of  a 
matched  filter-envelope  detector  combination  in  terms  of  the  ambiguity 
function  for  continuous  wave  (CW)  and  frequency  modulated  (FM)  input 
signals  when  these  signals  were  characterized  by  time  and  frequency 
spreading.  In  his  studies,  ambiguity  function  calculations  were  made  for 
both  noise  and  reverberation  backgrounds.  The  problem  of  correlation 
degradation  in  passive  systems  as  a  result  of  track  and  motional 
instabilities  in  a  transiting  source  (assumed  to  be  a  narrowband  signal 
source)  was  examined  by  Gerlach  [3].  In  his  work,  the  integration  time  to 
optimize  processor  gain  is  derived  for  signal  detection  with  an  incoherent 
noise  background.  Robertson  [4]  determined  the  performance 
characteristics  of  a  matched  filter-envelope  detector  for  multiple  CW 
pulses  in  narrowband  noise.  He  showed  the  improvement  in  receiver 
operating  characteristic  (ROC)  curves  as  several  independent  samples  of 
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the  envelope  of  the  combined  signal  and  noise  are  averaged  in  making  one 
attempt  at  detection.  Roberts  [5  ]  considered  the  detection  of  a  sinusoid 
with  signal  known  except  for  phase  which  was  described  by  a  unimodal 
distribution.  Kincaid  [6,7]  determined  optimum  waveforms  for  correlation 
detection  for  both  noise-limited  and  reverberation-limited  conditions. 

Most  of  the  significant  research  efforts  cited  here  were  based  on  the  more 
fundamental  radar  work  of  Marcum  [8],  Swerling  [9],  and  Schwartz  [loj. 

Each  considered  the  multiple  pulse  performance  of  matched  filter-envelope 
detectors  using  different  signal  fluctuation  models. 

In  this  report,  the  performance  of  the  quadrature  receiver  is 
theoretically  measured  in  terms  of  output  signal-to-noise  ratio  (SNR)  when 
the  input  signal  is  postulated  to  be  narrowband  with  Rayleigh-distributed 
amplitude  fluctuations  and  uniformly  distributed  phase  fluctuations  in 
time  and  the  background  noise  is  white  with  zero-mean.  As  the 
correlation  characteristics  of  the  input  signal  are  varied,  the 
degradation  in  the  receiver  output  SNR  is  shown  relative  to  the  case  when 
the  input  signal  is  coherent.  For  the  problem  considered  in  this  report, 
the  constant  signal  frequency,  f,  is  assumed  to  be  a  known  parameter. 
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SIGNAL  AND  NOISE  MODELS 

In  this  analysis,  the  received  waveform,  r(t),  in  Fig.  1  consists  of 

signal-plus-noise  or  noise  only.  The  noise,  n(t),  is  assumed  to  be 

white  with  zero-mean.  The  two-sided  noise  power  spectral  density  is 

N  /2  with  the  autocorrelation  function  equal  to  (N  /2)  5(t-t'); 
o  o 

hence,  the  noise  is  a  stationary  random  process.  The  signal,  s(t), 
postulated  to  be  narrowband  and  Gaussian,  has  the  following  mathematical 
representation: 

s(t)*V(t)cosCait+e(t)3  ,  for  0<t<T,  (1) 

where  V(t),  the  envelope,  is  a  Rayleigh  distributed  random  process  and 
0(t),  the  phase,  is  a  uniformly  distributed  random  process.  The  radian 
frequency,  oi,  is  assumed  to  be  known  and  constant  over  the  signal  duration, 
T.  Written  in  terms  of  in-phase  and  quadrature  components,  Eq. (1)  becomes 

s(t)*x(t)cosuit-y(t)  sinut,  for  0<t<T,  (2) 

where  x(t)«V(t)cos8 (t) ,  y(t)«V(t)sin8(t) .  x(t)  and  y(t)  are  the 
amplitudes  of  the  in-phase  and  quadrature  components,  respectively. 

In  underwater  acoustics  this  model  is  frequently  used  to  approximate 
the  signal  that  is  produced  when  many  overlapping  direct  path  and  surface 
path  returns  combine  after  being  reflected  from  a  distributed  highlight 
source.  Although  each  path  return  has  a  random  but  fixed  amplitude  and 
phase,  the  resultant  multipath- laden  signal  has  an  envelope  and  phase  that 
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are  time- varying  random  processes,  if  many  path  returns  combine  to 
produce  a  resultant  signal  like  the  one  expressed  in  Eq.(l),  the  Central 
Limit  Theorem  indicates  that  x(t)  and  y(t),  the  resultant  in-phase  and 
quadrature  components,  are  zero-mean,  Gaussian,  random  processes;  hence, 
the  envelope  V(t)  may  be  shown  to  be  a  Rayleigh  random  process  and  the 
phase  0(t),  a  uniform  random  process  [ll]. 


DERIVATION  OF  RECEIVER  OUTPUT  SNR 

In  this  section,  an  expression  is  derived  relating  the  output  SNR  of 
the  quadrature  receiver  at  time  T  to  the  input  waveform.  Beginning  with 
r (t)=s(t)+n(t)  in  Fig.  1,  it  is  evident  that  the  correlator  outputs  are 


■L 


q„(T)»  /  s(t)sin(ut)dt+ /  n(t)sin(wt)dt 


A”' 


(3) 


and 

T  T 

q_(T)*/'  s(t)cos(wt)dt  +  /  n(t)cos(u>t)dt, 

c  Jo  Jo 

where  E[qc(T)]  *E[qg(T)] -0.  (The  symbol  e[*]  indicates  expected  value.) 

It  is  assumed  for  simplicity  that  the  signal  and  noise  are  statistically 
Independent  stationary  random  processes.  Squaring  q  and  q  and  taking  the 
expected  value  of  their  sum,  Q(T),  yields 


T  T  T  T 

e[Q(T)] mj  J  Rg(t-t')cosw(t-t')dtdt '+J  j  Rn(t-t,)cosu>(t-t,)dtdt',  (4) 
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where  R  (t-t')»E[s(t)s(t')] ,  and  R  (t-t* ) »E[n (t)  n(t' )]  .  For  convenience, 
3  n 

R  (t-t*)  can  be  written  in  terras  of  the  autocorrelation  functions  of  x(t) 
s 

and  y(t)  by  using  s(t)«x(t)coswt-y  (t)  siruut,  and  knowing  that  for 

the  narrowband  representation  in  Eg. (1)  R  (t-t')*R  (t-t*)  and 

x  y 

R  (t-t')«  -R  (t-t')  as  shown  in  [l2],  Hence,  after  making  these 
xy  yx 

substitutions  and  applying  trigonometric  identities  the  autocorrelation 
function  of  s(t)  is  written  as 


R  (t-t')-R  (t-t')cosu(t-t*)+R  (t-t'UinaKt-t') .  (5) 

s  x  xy 

Substituting  Eq.{5)  into  Eq. (4)  and  ag.\in  applying  trigonometric  identities 
gives 


.  T  T 

E[Q(T)]  »  *u  y* ^dtdt  + 


.  T  T 

tf.fr  (t-t'  )cos2u)(t-t' )  dtdt' 


,  T  1 

+  Vo/, 


R  (t-t' )sin2u(t-t') dtdt 
xy 


T  ’I 

'I/, 


R  (t-t'  )cosaj(t-t'  )dtdt' 
n 


(6) 


Following  Turin’s  development  of  the  SNR  for  a  matched  filter-envelope 
detector  matched  to  a  non-coherent  signal  in  noise  [l3],  the  double  frequency 
terms  in  Eq.(6)  are  neglected.  In  fact,  when  the  product  of  T  and  u>  is  large, 
as  it  often  is  in  sonar  applications  (in  the  kilohertz  range),  these  double 
frequency  terms  vanish  (see  appendix).  Using  the  noise  model  postulated  in 
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the  previous  section  to  evaluate  the  integral  in  Eq.(6)  with  Rn(t-t')» 

Eq. (6)  becomes 

T  T  T  NT 

E[Q(T)]s|y  j  Rx(t-f)dtdt,+ (7) 

Since  e[q(T)]  represents  the  average  receiver  output  power  at  t*T,  the 
approximate  output  SNR  at  t*T  is 

7  T  T 

SNRS_±_/  /  R  (t-f)dtdf.  (8) 

o  Jo  Jo 

Papoulis  [14]  shows  how  this  double  integral  can  be  simplified  to  the 
following  single  integral 

SNRSN_|  (l"  T  )Rx(T)dT-  (9) 

This  result  explicitly  relates  the  output  SNR  to  the  noise  power  spectral 
density,  Nq/2;  the  integration  time,  T;  and  the  in-phase  component  of  the 
signal  (or,  equivalently,  the  quadrature  component). 
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SNR  VS.  SIGNAL  CORRELATION 


The  output  SNR  will  now  be  evaluated  £or  signals  with  different  correla¬ 
tion  characteristics.  First,  consider  the  case  where  the  amplitude  and 
phase  of  the  signal,  s(t),  are  constant  and  unknown  over  the  duration,  T 
(i.e.,  V  and  6  are  constant).  The  result  for  this  coherent  signal  case  is 
well  known  and  is  presented  by  Whalen  Cl]  in  his  treatment  of  signals  with 

random  phase  and  amplitude.  In  this  case,  the  autocorrelation  function 
2  2  2 

R  (T)  becomes  E(V  cos  6)  or  Aq/  where  Aq  is  defined  as  the  Rayleigh 
parameter  describing  the  Rayleigh  variable  V.  Hence,  Eg. (9)  reduces  to 


SNRJ»~  .  (10) 

o 

This  result  agrees  with  Whalen's  development. 

Next,  x(t)  is  considered  to  be  a  Gauss-Markov  process  where  the 
autocorrelation  function  is  Rx(T)»A^e  a^Tlf  where  1/a  is  defined 
as  the  correlation  time.  Substituting  into  Eq.(9)  and  performing  the 
integration  yields 


— aT 

2(e  -1+aT) 

A  1 

No 

(aT)2 

(11) 
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From  Eq.(ll),  as  a  ♦  ®,  or  as  the  correlation  of  x(t)  with  itself 

2 

diminishes,  the  SNR  -*■  0.  Also,  as  a  0,  the  SNR  *  A  T/N  ;  this  result 

o  o 

agrees  with  the  coherent  signal  case  considered  earlier. 

Finally,  the  case  where  x(t)  is  a  Gauss-peaked  process  with 
Rx(T)*A^e”a^  coso>o|t|  is  considered.  Here  Eq.  (9)  becomes 


where  c*a+ juQ ,  a  complex  variable.  When  uo*0,  c  equals  the  real  value 

a;  x(t)  becomes  a  Gauss-Markov  process;  and  Eq. (12)  reduces  to  Eq. (11) . 

Since  Eq. (12)  is  in  effect  a  composite  of  each  of  the  other  cases,  it  is 

presented  graphically  in  Fig.  3.  In  the  figure,  the  SNR  is  normalized  with 
2 

respect  to  AqT/Nq ,  the  result  obtained  for  a  coherent  input  signal,  and 
plotted  in  decibels  as  a  function  of  aT  with  u>oT  a  parameter. 
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RESULTS 

Figure  3  describes  the  performance  of  the  quadrature  receiver  as 
certain  signal  parameters  are  varied.  It  should  be  remembered  that  the 
integration  time,  T,  remains  fixed  and  equal  to  the  signal  duration;  hence 
the  only  parameters  that  are  varied  are  the  correlation  time,  1/a,  and  the 
radian  frequency  of  oscillation,  u)Q,  of  the  autocorrelation  functions 
describing  x(t)  and  y(t). 


In  figure  3,  the  oj  T»0  curve  demonstrates  the  receiver's 

o 

performance  when  x(t)  and  y(t)  are  Gauss-Mar kov  processes;  the  origin 

aT 


Figure  3.  Normalized  Output  SNR  for  Quadrature  Receiver 
When  x(t)  and  y(t)  Are  Gauss-peaked  Processes 
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corresponds  to  the  coherent  input  signal  case.  The  other  curves  represent 

the  case  where  x(t)  and  y(t)  are  Gauss-peaked  (i.e.,  exponential-cosine) 

randan  processes.  It  is  interesting  to  observe  the  behavior  of  the  slope 

of  the  curves  at  aT*0.  The  u  T*0  and  n>  T»3  curves  have  a  negative 

o  o 

slope  at  this  point  and  are  nonotonic  decreasing  functions  of  aT.  The 
(uqT«4  and  o>oT«5  curves  have  a  positive  slope  at  the  origin  and 
peak  before  again  becoming  monotonic  decreasing  functions.  It  can  be 
shown  that  the  turnover  point  for  the  slope  at  aT*0  occurs  woT*it ; 
here,  the  slope  is  zero. 

Finally,  it  is  worthwhile  to  note  that  there  appears  to  be  a  value  of 
aT  beyond  which  the  SNR  is  independent  of  (t>oT. 
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SUMMARY 

In  active  sonar  systems,  when  the  transmit  signal  is  a  narrowband  CW 
pulse,  a  combination  of  non- resolvable  multipath  returns  from  a  reflector 
degrade  the  detection  performance  of  a  quadrature  receiver.  This  well- 
known  receiver  is  optimum  when  the  input  waveform  consists  of  a  coherent 
signal  with  a  white  noise  background.  This  combination  of  non-resolvable. 
multipath  returns  causes  the  amplitude  and  phase  of  the  signal  to  become 
random  functions  of  time.  As  a  result,  the  signal  diminishes  in  coherence 
and  the  receiver's  performance  is  degraded.  The  work  presented  in  this 
report  has  theoretically  measured  the  degradation  in  a  quadrature 
receiver's  output  SNR  using  a  narrowband  signal  model  with 
Rayleigh-distributed  envelope  fluctuations  and  uniformly  distributed  phase 
fluctuations  over  the  signal  duration.  Several  variations  in  the 
autocorrelation  function  of  the  signal  were  considered.  After  deriving  an 
expression  for  the  output  SNR  of  a  quadrature  receiver 
showing  its  dependence  on  the  autocorrelation  function,  Rx(t),  of  the 
in-phase  (or  quadrature)  amplitude,  the  output  SNR  was  determined  for 
several  variations  of  Rx(t).  For  comparison  purposes,  the  SNR  was 
normalized  to  that  obtained  when  the  input  signal  is  coherent. 
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APPENDIX 


DISCUSSION  OF  DOUBLE  FREQUENCY  TERMS 


This  development  demonstrates  why  the  double  frequency  terms  in 
Eq. (6)  can  be  neglected.  To  begin,  it  is  known  from  the  Wiener-Kinchine 
theorem  that 


where  S^Coj)  is  the  power  spectral  density  of  the  random  process  x(t). 
Hence, 
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Evaluating  the  inner  integrals  yields 
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(A3) 


Recognizing  that  Eq. (A3)  can  be  written  in  terms  of  sine  functions  gives 
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It  is  evident  that  Eq. (A4)  is  small  when  Tea.  >>  1. 
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